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ABSTRACT 

Context. Synthetic model atmosphere calculations are still the most commonly used tool when determining precise stellar parame¬ 
ters and stellar chemical compositions. Besides three-dimensional models that consistently solve for hydrodynamic processes, one¬ 
dimensional models that use an approximation for convective energy transport play the major role. 

Aims. We use modern Balmer-line formation theory as well as spectral energy distribution (SED) measurements for the Sun and 
Procyon to calibrate the model parameter a that describes the efficiency of convection in our ID models. Convection was calibrated 
over a significant range in parameter space, reaching from F-K along the main sequence and sampling the turnoff and giant branch 
over a wide range of metallicities. This calibration was compared to theoretical evaluations and allowed an accurate modeling of 
stellar atmospheres. 

Methods. We used Balmer-line fitting and SED fits to determine the convective efficiency parameter a. Both methods are sensitive to 
the structure and temperature stratification of the deeper photosphere. 

Results. While SED fits do not allow a precise determination of the convective parameter for the Sun and Procyon, they both favor 
values significantly higher than 1.0. Balmer-line fitting, which we find to be more sensitive, suggests that the convective efficiency 
parameter a is » 2.0 for the main sequence and quickly decreases to « 1.0 for evolved stars. These results are highly consistent with 
predictions from 3D models. While the values on the main sequence fit predictions very well, measurements suggest that the decrease 
of convective efficiency as stars evolve to the giant branch is more dramatic than predicted by models. 

Key words. Stars: atmospheres - Stars: fundamental parameters - Stars: late-type - Stars: Balmer lines 


1. Introduction 


One-dimensional (ID) stellar model atmospheres are widely 
used for determining basic physical parameters, such as effective 
temperatures, surface gravities, luminosities, masses, and chem¬ 
ical compositions of stars at different evolutionary stages. They 
also serve as boundary conditions for stellar evolution codes. 

While codes for 3D model atmospheres self-consistently 
solve for convection in the atmospheric layers© D codes apply 
a simplified appr oach such as described bvlBohm-Vite nse 1958 
(BV hereafter) or lCanuto & Mazzitellill 199 ill 19921 (CMA here- 
after). As these formulations do not self-consistently determine 
the convective flux, at least one calibration parameter is required: 
the a-parameter, which sets the efficiency of convection. 

Recent calibrations of the ID a-parameter based on 3D the¬ 
oretical model atmospheres have been presented by Magi c et alJ 
©0141 (hereafter. Magic 14) and earlier bv iLudwig et all ( 19991) . 
The primary goal of this paper is to compare observational 
constraints on acMA , the a -parameter in the formulation of 
ICanuto & Mazzitelh] d 1 9921) . to theoretical predictions of 3D 
models. This will allow us to use a proper ctcma for ID model 
atmospheres across a wide range of stellar parameter space. 


Send offprint requests to: F. Grupp; e-mail: frank@grupp-astro.de 


We use the approach as iFuhrmann et al.lfl993l (FAG93 here¬ 
after), using the dependence of Balmer line profiles on the at¬ 
mospheric structure, which we call the ffcMA-parameter. For this 
test we use a sample of stars (including the Sun) with stellar 
parameter start values determined independently of model at¬ 
mospheres, for example, by astrometrical methods. As another 
tool for determining acMA , we use stellar spectral energy dis¬ 
tributions (SED) in the visible and near-UV spectral range. This 
method is of course limited to objects with good absolute cali¬ 
bration of the SED. We here use the Sun and Procyon. 

On the model side, the adopted model atmosphere code is the 
plane-parallel, chemically homogeneous, local thermodynami¬ 
cal equilibrium (LTE) code MAFAGS-OS presented by Grunnl 
( 2004a}) . Stellar convection is treated in MAFAGS-OS accord¬ 
ing to the formalism of CMA. The authors take into account a 
spectrum of turbulent eddies with different sizes, which is more 
advanced than the one single eddy in the mixing-length the¬ 
ory (MLT; iBohm-Vitensell 195 8). The molecular weight within 
the convective zone is variable, a detail neglected by MLT. 
Nevertheless, it is a parametric and simplified approach, not a 
self-consistent numerical formalism. MLT furthermore assumes 
that the turbulence is incompressible and sets the mixing-length 
A = o7/ p , where H p denotes the local pressure scale height, 
and leaves the so-called mixing-length parameter a (hereafter. 
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otvilt) as a free parameter. Therefore, calibrations by comparing 
the models with the observational data are necessary. 

Some previous works suggest that ctmlt is greater than unity 
for the Sun and var i es wit h stellar parameters. For example, 
iMiglio & MontalMnl ( 2005 ) constrained »mlt of the a Cen- 
tauri A+B systems using a method from asteroseismology and 
found an otmli for component B higher by 10% than that of A. 
Moreover, ctmlt was_ found to be strongly dependent on stel¬ 


lar mass ( e.. 

[Ferraro et al 


g.. [Yildiz et al.l 12006 ), but not on metallicity (e.g., 

_ _7l2006l). In contrast, us ing the solar-like stars in 

the Kepler field. iBonaca et al.1 (1 201 2l) indicated that ctmli corre¬ 
lates more significantly with metallicity than 7’ e ii or logg. With 
2D hydrodynamical models of dwarf stars. iLudwig et alJ(il999l) 
showed that ctmlt decreases with increasing T e g and log #, and 
this was confirmed by 3D simulations of lTramnedachl ( 2007 ). 

ICanuto & Mazziteifil ( 1992 ) introduced A = az by taking 
into account local effects, where z is the distance to the top of the 
convective zone, and its value is clo se to the poly trope H„. a = 1 
correspond to the nonlocal A = z in lCanuto & Mazzitelli ( 199 D. 
In previous versions of the MAFAGS-OS program, the a accord¬ 
ing to the CMA theory (hereafter, ctcma) is fixed to 0.82, because 
the evolutionary stage of the Sun, as determined by its internal 
structure, can be well reproduced by this value (iBernkonfj 1998l) . 
However, this is only one calibration point in parameter space, 
there is no evidence that stars throughout the whole H-R diagram 
necessarily have the same »cma values. 

In this paper, we analyze a sample of extensively studied 
"standard stars” Q with different chemical compositions and evo¬ 
lutionary stages. All of them have angular diameter data ob¬ 
tained with modern interferometers, which - together with dis¬ 
tance measurements from Hipparcos - enable direct measure¬ 
ments of stellar radii, effective temperature, and surface gravi¬ 
ties. The variation of ctcma across the H-R diagram can therefore 
be constrained. 


2. Stellar samples and basic stellar parameters 

Our sample consists of ten dwarf and turnoff stars as well as 
one giant, with effective temperatures ranging from 4600 K to 
6600 K and metallicities ([Fe/H]) from -2.60 to +0.60. All 
objects have been observed with the Fiber Optics Cassegrain 
Echelle Spectrograph (FOCES: lPfeiffer et al,lll998h on the 2.2m 
telescope at Calar Alto Observatory. The optical spectra have 
signal-to-noise ratios of >200 and resolving powers (R = A/ Ad) 
of «60,000. Data reduction was_performed using the FOCES 
EDRS software package (see lPfeiffer et al.ll998l) . Normalization 
was performed manually. The sample is shown in the physical 
Hertzsprung-Russell-Diagram(HRD) in Fig.Q] Interferometric 
diameter determination for eight of the program stars is avail¬ 
able (see Table[T|). 

To minimize the different instrumental response effects, the 
spectrosco pic p arameters (7/n, logg, [Fe/H] and / : j were taken 
from [Fuhrmannl (fl998l 120041) . who derived them with the same 
FOCES spectra as we used in this study. The available in¬ 
terferometric observations provide accurate bolometric fluxes 
CTbql ) with typical uncertainties of ~2% (e.g., iBovaiian et alJ 
120121) . Combined with the interferometrically measured stellar 
angular diameter (0) measurements with typical relative uncer¬ 
tainties of 1~3% for our samples, the uncertainties of T e g = 


1 FG would like to suggest defining a "standard star" as a star that - 

as soon as it is studied in detail - is found to be pretty "non-standard". 

Procyon and HD 140283 are stars that fit this category very well. 
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Fig. 1. HRD of the program stars. Color encodes metallicity, size en¬ 
codes stellar diameter measurements. 


2341((Fbol/# 2 * ) 1 ^ 4 are lower than 1.5%, which is precise enough 
as cross-validations of the spectroscopic data. 

Therefore, whenever possible, we used interferometric 
and 0 given by recent measurements with VLTI or CHARA to 
compose interferometric parameters. The spectroscopic and in¬ 
terferometric parameters are summarized in Table Q] 

The log g of stars that have interferometric measurements 
parameters was derived using the direct angular diameter of 
the star, its distance from the Hipparcos parallax (Ivan Leeuwen! 
120071) . and its mass, either from the astrometric measurements 
of a binary companion (for Procyon A), or by fitting the stellar 
evolutionary track in the H-R diagram. The parameters metallic¬ 
ity ([Fe/H]) and micro-turbulence velocity (£■'), which cannot be 
obtained by interferometric means, were kept the same as those 
in the spectroscopic group. The solar abundance mixture was 
taken from iLodders et all ( 2009 ). For the five low-metallicity 
stars in our sample with [Fe/H] < -0.60 (HD 6582, [Fe/H] 
= -0.83; HD 103095, [Fe/H] = -1.35; HD 122563, [Fe/H] = 
-2.60; HD 45282, [Fe/H] = -1.50; and HD 140283, [Fe/H] 
= -2.38), Q'-elements enhancements of [a/Fe] = +0.4 were 
adopted. 

A model of Vega was calculated to consistently derive the 
theoretical zero-point of the spectrophotometric data. 


3. Model atmospheres 

For each set of stellar parameters, model atmospheres were cal¬ 
culated! usingjhemodelattnosphere code MAFAGS-OS (Grupp 
2004 albl ; iGrupp et alJl2009l) . The opacity-sampling (OS) code 
is based on the opacity distribution function (P DF) version o f 
T. lGehrenl (i 1 9791) . In the reprogrammed version nf lReileld 19871) . 
the code relies on the following basic assumptions: 

- Plane-parallel ID geometry. 

- Chemical homogeneity throughout the atmosphere. 

- Hydrostatic equilibrium. 

- Convection is treated accordin g to the formalism of 
ICanuto & Mazzitellild 199 111 19921) . 

- No chromosphere or corona. 

- Local thermodynamical equilibrium. 

- Flux conservation throughout all 80 layers. 

While these assumptions can break down for hot stars, stars with 
very extended atmospheres, and the coolest stars, they are valid 
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Table 1 . Best-fitting qtcma 


HD 

feff(K) 

@L D 

Ref. 

logg“ 

[Fe/H]“ 


The best 

H a Up 

Q'CM 

H r 

H ff 

*; e / d - af ) 

HjS 

H y 

Interferometric and IRFM 

Sun 

5780±40 

- 


4.44 

0.00 

5780 

2.0 

>1.0 

>1.0 

6.05 (244) 

- 

- 

Procyon 

6530±50 

5.448±0.053 

Kervella et al. (2004) 

3.98 

-0.05 

6524 

1.9 

>1.0 

>1.0 

0.94 (309) 

- 

- 

10700 

5380±50 

2.078±0.031 

Di Folco et al. (2004) 

4.53 

-0.49 

5330 

1.8 

>1.0 

>1.0 

0.72(195) 

- 

- 

103095 

4820±100 

0.679±0.015 

Creevey et al. (2012) 

4.59 

-1.30 

5100 

2.4 

1.9 

>0.5 

4.22 (648) 

1.84 (217) 


39587 

5960±50 

1.051 ±0.009 

Bovaiian et al. (2012) 

4.48 

-0.07 

5960 

1.9 

>1.0 

>0.5 

1.11 (187) 

- 

- 

6582 

5260±85 

0.972±0.009 

Bovaiian et al. (2012) 

4.51 

-0.83 

5370 

1.7 

>0.5 

- 

0.58 (112) 

- 

- 

217014 

5800±55 

0.748±0.027 

Baines et al. (2008) 

4.28 

+0.20 

5780 

2.0 

>1.0 

>1.0 

1.08 (236) 

- 

- 

22049 

5120±75 

2.148±0.029 

Di Folco et al. (2004) 

4.61 

-0.09 

5045 

2.2 

>1.0 

- 

0.83 (252) 

- 

- 

140283 

5780±50 

- 

Casagrande et al. (2010) 

3.70 

-2.38 

5780 

1.7 

1.0 

1.0 

1.00(201) 

0.60 (547) 

1.13 (264) 

122563 

4600±50 

0.948±0.012 

Creevey et al. (2012) 

1.60 

-2.60 

4600 

1.0 

0.5 

0.5 

0.46 (533) 

2.50 (236) 

2.43 (80) 

45282 

5300±50 

- 

Casagrande et al. (2010) 

3.07 

-1.50 

5310 

1.0 

0.5 

- 

0.96 (574) 

0.71 (124) 

- 

Spectroscopic 

Procyon 

6530 

- 

This work 

3.96 

-0.05 

6530 

1.9 

>1.0 

>1.0 

- 

0.99 (309) 

- 

103095 

5110 

- 

Fuhrmann (1998) 

4.66 

-1.35 

5110 

2.4 

1.9 

>0.5 

- 

1.05 (648) 

1.07 (217) 

39587 

5920 

- 

Fuhrmann (2004) 

4.39 

-0.07 

5920 

1.5 

>0.5 

>0.5 

- 

0.96(187) 

- 

6582 

5390 

- 

Fuhrmann (2004) 

4.45 

-0.83 

5390 

2.3 

>1.0 

- 

0.72(112) 

- 

- 

217014 

5790 

- 

Fuhrmann (1998) 

4.33 

+0.20 

5790 

2.0 

>1.0 

>1.0 

0.93 (236) 

- 

- 

22049 

5050 

- 

Fuhrmann (2004) 

4.54 

-0.09 

5050 

2.2 

>1.0 

- 

0.80 (252) 

- 

- 


Notes. “ References of the adopted logg and [Fe/H] are in the text. T* s denotes the best effective 
temperatures obtained from the fitting procedure. 


for temperatures from 4000 through 15000 K and for gravities 
from the main sequence down to log g « 0. 

The model was iterated until the flux was constant (< 0.5% 
level) throughout the atmosphere. 

3.1. Convection and atmospheric structure 

Figure[2]shows the influence of the convective efficiency param¬ 
eter ffcMA on temperature stratification, gas pressure, and elec¬ 
tron pressure of a solar atmospheric model. While a low value 
of a cma results in low convective flux, high »cma leads to an 
increase of the convective flux and - as the total flux is preserved 
through all layers - lower radiative flux in the convection zone 
of the stellar atmosphere. Low radiative flux leads to lower tem¬ 
perature gradients. Therefore, the high-crcMA case shows a lower 
temperature in the inner region of the stellar atmosphere. 

For our program stars, increasing q-cma from 0.5 to 2.5 leads 
to decreasing temperatures in the inner atmosphere by up to 

1000-1300 K at log(r 50 oo) - 1.5. 

The effect of a variation of ctcma on the solar SED is shown 
in Fig. [3 Flux changes of up to 5% of the solar SED in the UV 
and of 1% in the near-IR region are notable. There are some 
node points at 550, 1350, and 2050 nm around which the fluxes 
are insensitive to mixing-length parameters, while the continuum 
of the SED between each two shows opposite tendencies with 

O'CMA- 

The trends in model structures and fluxes are similar for other 
stars in our sample, with the only exception of the A0V star 
Vega. As Vega has no convection zone reaching the stellar at¬ 
mosphere, no change with ac ma is present in the stellar flux. 


4. Spectral energy distributions 

In the following section we directly compare the computed ab¬ 
solute fluxes to the observational data to calibrate acMA- Abso¬ 
lute flux data are only available for a few stars. For our sample, 
reliable absolute SED data are only available for the Sun and 
Procyon. 

4.1. Near-UV - visible - near-IR domain 

iNeckel & Labs 1 1984 1 measured solar irradiance between 3300 
and 12500 A by combining the absolute fluxes of high-resolution 
spectra obtained with the Fourier Transform Spectrometer (FTS) 
at the McMath Solar Telescope on Kitt Peak. The spectra were 
averaged every 10 A between 3305 and 6295 A, and every 20 A 
between 6310 and 8690 A. 

iThuillier et al.l (2003) presented the solar irradiance from 
2000 A to 24000 A obtained with the SOLSPEC and SOSP spec¬ 
trometers mounted on space missions. However, for wavelength 
ranges above 8700 A, their data are sampled with a rate different 
from that in visible and UV domains. Therefore we only adopted 
the spectra below 8700 A to compare with our solar fluxes. 

Figure® sh ows the data and absolut e deviations of sola r ir¬ 
radiance o nNeckel & Labsl (i 19841) and lThuillier et al.l (120031) to¬ 
gether with the computed solar fluxes for three different »cma 
values of 0.5, 1.5, and 2.5. A visual check prefers a fit between 
blue (a C MA = 2.5) and green (a C MA = L5). 

The left and right panels of Fig.[5]compare the measured and 
observed fluxes starting from 4500 A and 5500 A. The strong de¬ 
pendency on the blue end of the comparison window does not 
allow precisely determining a cm a from these data. Neverthe¬ 
less, data points with values below 1.0 can be excluded from the 
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Fig. 4. Comparison of solar absolute fluxes { solid dots ) to synthetic flux levels with o cma = 0.5 {red), 1.5 {green), and 2.5 (blue). 

Observational data were taken from lNeckel & Labsld 19841) {upperpanel) and lThuillier et al J d2003l) {lower panel). The absolute deviation (F mode i - 
F 0 b s ) is plotted in the lower part of each panel. Offsets of 0.05 were added to improve the visualization. 


Neckel data. The data of lThuillier et al.l (120031) would only allow 
for very low convective efficiency if the blue region is ignored. In 
general, we find an overestimation of the solar flux in the regime 
below 5000 A. This can be caused by missing atomic or molec¬ 
ular opacity. Variations of convective efficiency do not change 
this region significantly, thus can be excluded as the cause of 
this shortcoming. 

4.2. Ultraviolet domain 

The ultraviolet spectrum of ProcyonA has been obtained with 
the Space Telescope Imaging Spectrograph (STIS) onboard the 
Hubble Space Telescope (HST) with high-quality calibration in 


the range from 220 nm up to 405 nm. As the STIS spectra span 
several echelle orders and overlap between each pair of orders, 
we averaged bins of 2 nm width for the STIS measured data and 
MAFAGS-OS model calculations to allow a most direct com¬ 
parison. This comparison is shown in Fig. 0 for a C MA rang¬ 
ing from 0.5 to 2.5 for a model with a 7’ e ff/logg/[Fe/H]/f = 
6530/3.98/-0.05/2.10. MAFAGS-OS underestimates the fluxes 
by at most ~35% for Mg n H & K at 280 nm. The average devi¬ 
ation is within 4%. 

Figure |T| displays the 2D plots of x 2 versus »cma (x axes) 
and the scaling factor of computed flux (y-axes) for four models 
with different temperatures. The difference of ~ 70 K in F e g- 
causes a difference of ~4% in the total fluxes. The scaling factors 


Article number, page 4 ofl 141 





























Wu et al.: CM A calibration 




0.0075 

0.0074 

0.0073 

0.0072 

0.0071 

0.0070 

0.0069 

0.006g 




Fig. 5. Relations of^ 2 of observational and computed fluxes versus q-cma for the solar irradiance of iNeckel & Labs! ( t 1984l f and iThuillier et al] 
( t2003f) . Observational data between 4500 and 8700 A (left plots) and between 5500 and 8700 A (right plots) were taken into account. 


were introduced to account for the offsets between the computed 
fluxes and the observational data to first-order approximation. 
The relative uncertainty (~ 1%) of the angular diameter can lead 
to an offset of up to 2% of the observational SED. Therefore the 
±3% of the scaling factor on the total flux is reasonable. The 
derived «cma for each model was estimated by minimizing the 
X 2 with the observational SED over the 2D plane. We found that 
model (c) and (d) led to smaller^ 2 with the observational data, 
which favors higher T e g for Procyon A. Models (b), (c), and (d) 
led to quite consistent ctcma — 1.4 to minimize ^ 2 , while model 
(a) favors smaller a cma — 1.3. 


5. NLTE calculations of Balmer-line profiles 

Our analysis is based on the nonlocal thermodynamical equilib- 
ri um (NLTE) line formation for H i using the method described 
in iMashonkina et al.l (120081) . In brief, the model atom includes 
levels with principal quantum numbers up to n < 19. We checked 
the influence of inelastic collisions with neutral hydrogen atoms, 
as computed following ISteenbock & Holwegeri d!984l) . on the 
statistical equilibrium (SE) of Hi and resulting profiles of the 
Balmer lines in the models. The differences in normalized fluxes 
between including and neglecting H+H collisions are smaller 
than 0.0005 of the continuum flux for the core-to-wing transi¬ 
tion region that is most sensitive to T e g variations. Therefore all 
the NLTE calculations were performed with pure electronic col¬ 
lisions with the exception of a very metal poor (VMP) giant at¬ 


mosphere 4600/1.60/—2.6 and a metal-poor subgiant atmosphere 
5777/3.70/-2.38, for which S#=1.0 and 2.0 were adopted. 


The coupled radiative transfer and SE eq uatio ns were solved 
with a revised version of the code DE TAIL ( Butlei^&Giddingsl 
119851) . The update was described in IMashonkina et alJ (12008). 
The obtained departure coefficients were then used by the code 
SIU ( Reel/ 199 1) to calculate the synthetic line profiles. The 
absorption profiles of Balmer lines were computed by convolv¬ 
ing the profiles resulting from the thermal, natural, and Stark 
broadening, as well as self-broadening. The Stark profile was 
based on the unified theory as de veloped by Vidal et akl (I1970L 
1973, VCS). The calculations of Stehlel (11994 based on a dif¬ 
ferent theory agreed reasonably well with the data of VCS. 
For self-broadening, we applied the Lorentz profile with a half- 
width computed using the cross-section and velocity paramete r 
from lBarklem et akl (120001) . As shown bv lBarklem et al.l ((2000), 
"the complete profiles obtained from overlapping line theory are 
very closely approximated by the p-d component of the relevant 
Balmer line". This gives ground to apply the impact approxima¬ 
tion to describe the self-broadening effect on H,,, H^, and H,,. 


We selected the model 6530/3.96/-0.05 to show the behavior 
of the departure coefficients bj = rtf [AV ‘/n\ A H of the relevant lev¬ 
els of Hi (Fig. [8]). Here, «^ LTE and n ETE are the SE and thermal 
(Saha-Boltzmann) number densities, respectively. Very similar 
behavior of bj was found for all the model atmospheres investi¬ 
gated in this study. The ground state and the n — 2 level keep 
their thermodynamic equilibrium level populations throughout 
the atmosphere. Departures from LTE for the n = 3 level are con- 
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Fig. 2. Temperature (T), pressure ( P ), and electron pressure (PA versus 
optical depth log(r 50 oo) with three different mixing-length parameters 
ctCMA in the MAFAGS-OS models of Sun. 


trolled by H„. In the layers where the continuum optical depth 
drops below unity, H ff serves as the pumping transition, which 
results in an overpopulation of the upper level. For instance, 
b(n = 3) =* 1.017 around logrsooo = -0.57, where the H ff core- 
to-wing transition is formed. Starting from log T 5000 — -2.8 and 
farther out, photons that escape from the Hq, core cause an un¬ 
derpopulation of the n = 3 level. For H„, NLTE leads to a weak¬ 
ening of the core-to-wing transition compared to the LTE case 
(Fig.®. 

To compare them with observations, the computed syn¬ 
thetic profiles were convolved with a profile that combines in¬ 
strumental broadening with a Gaussian profile of 3.2 km s 1 to 



Fig. 3. Upper panel: Computed solar flux F A with different ffcMA- The 
spectrum is sampled every 50 A along the wavelength. The transmission 
curves of U, B, V, R. I in Bessell photometric system, and J, H, K s in 
2MASS system are also overplotted (from left to right). The residuals 
to fluxes with »cma = 0.82 are plotted in the lower panel. 



Wavelength (nm) 

Fig. 6. Comparison of absolute fluxes of Procyon A (solid lines) ob¬ 
tained with HST/STIS with MAFAGS-OS model fluxes with «cma = 
0.5 {red), 1.5 (green), and 2.5 (blue) in model (a). The fluxes were 
summed every 2 nm. The absolute deviation (F mode i - F obs ) were plotted 
in the lower panel. Offsets of 1.0 in the upper panel and 0.5 in the lower 
panel were added for clarity. 

4.6km s' for different stars and broadening caused by macro¬ 
turbulence with a radial-tangential profile of 4kms _1 . Rota¬ 
tional broadening was only taken int o ac count for Procyon, 
with Vsini = 2.6kms' 1 (lFuhrmannlll998l) . All the remaining 
dwarf and subgiant stars were assumed to be slow rotators, with 
V sin i < 1 km s' 1 . It is worth noting that the thermal broaden¬ 
ing of Balmer lines, with the most probable velocity of about 
V, = 10kms' 1 in the atmospheres of the investigated stars, 
is more pronounced than the broadening caused by rotation, 
macroturbulence, and instrumental profile together. Therefore, 
the uncertainty in the external broadening parameters only has 
an minor effect on the computed Balmer-line wings. 
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Fig. 7. X 2 versus q;cma and the scaling factor for the HST/STIS 
data and the computed absolute fluxes between 220 and 405 nm 
with models of different T t g for ProcyonA. From (a) to (d): r e ft = 
6520/6530/6550/6590 K. 



T 5000 


Fig. 8. Departure coefficients log bj for the five lowest levels of Hi as a 
function of continuum optical depth r 50 oo referring to A = 5000A in the 
model atmosphere 6530/3.96/-0.05. Tick marks indicate the locations 
of core-to-wing transition (0.85 to 0.95 in normalized flux) formation 
depths for H ff and Hg. 

6. Fitting procedure of Balmer-line wings 

Our analysis of Hq, is based on profile fitting in the wavelength 
range from 6473 to 6653 A for Procyon and from 6520 to 6608 
A for the remaining stars. We analyzed Hg in the wavelength 
range 4815 - 4908 A only for the metal-poor stars HD 103095, 
HD 122563, HD 45282, and HD 140283. For HD 122563 and 
HD 140283 we also analyzed H y profiles in the wavelength range 
from 4315 to 4365 A . Before fitting, we made the following 
preparations of the observed profiles for all stars. 

We excluded the Balmer line cores up to 0.9 of the rela¬ 
tive flux (Fj/Ff), because the core half-width of the hydrogen 
line cannot be described in cool stars with classical model atmo¬ 
spheres. Thus only the pressure-broadened wings were analyzed. 
As our aim is to fit the profiles of Balmer lines, all blended lines 



Wavelength (A) 


Fig. 9. Synthetic NLTE (continuous curve) and LTE (dashed curve) 
flux profile of H, computed with the MAFAGS-OS 6530/3.96/-0.05 
(ctcma = 2.3) model compared to the observed spectrum of Procyon 
(bold dots). 


should be removed from the observed profiles. The high resolu¬ 
tion of the employed spectra allows us to assume that unblended 
spectral windows do exist and can be used for the fitting. We 
identified spectral windows across the line that are expected to be 
free from blends in the stars of interest. This was done by means 
of calculating Balmer line profiles with and without metal lines. 
Then we created a mask that was free from any blends. Based 
on the examination of the spectra, this is reasonable at H 0 for 
all the reference stars, and at H /; for HD 122563, HD 45282, and 
HD 140283, and at H y only for HD 122563 and HD 140283. It is 
impracticable to define spectral windows without any blends at 
Hj S and H y in solar metallicity stars. 

The mask consists of a spectral window of pure H,, , H /; , and 
H y (when available) lines with noise alone. The same mask was 
applied for all stars, but minor modifications were made for each 
star individually. Modifications of the mask were necessary to 
reject stars in some specific cases, for example, if there was an 
ingress of atmospheric lines in the spectral window or for Pro¬ 
cyon, which has much stronger macroturbulence and V sin / than 
other stars. It is very important to clearly distinguish between 
noise fluctuations and weak metal lines in the spectra. In this 
context, we relied upon the spectral line database and inspected 
all Balmer line profiles. 

We also investigated the effect of spectral resolution on the 
choice of which spectral windows to mask. It is necessary to 
have a proper number of spectral windows between blended lines 
because only by using such a mask are we able to derive the same 
result regardless of whether we use high-resolution spectra or the 
spectra degraded to the resolution of our observations. This mask 
is quite reliable in the case of H„ for all stars, but for Hg and H y 
we were unable to define such a mask for all stars. This makes 
the results obtained from Hg and H y profiles less reliable. 

We also took into account the effect of rotational broaden¬ 
ing and broadening caused by macroturbulence on the choice 
of spectral windows of mask. While all our sample stars are 
slow rotators with V sin i < 1 km s 1 except for Procyon with 
V sin/ = 2.6 km s^ 1 , we were able to find a proper number of 
spectral windows between blending lines around H,, profiles. In 
the case of Hg and H y profiles, this effect becomes important 
because there are many blending lines in the two profiles. This 
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prevented us from defining a reliable mask for all stars, except 
for some metal-poor stars (HD 103095, HD 122563, HD 45282, 
and HD 140283), which only contain few blending lines in their 
Hyj and H y profiles. 

In this work we used a method based on a reduced^ 2 statis¬ 
tics. Statistical measurement of the goodness-of-fit was per¬ 
formed by comparing the theoretical profile and an observed 
profile through a x 2 minimization procedure (INousek & Shuei 
Il989h . The estimator of the fit goodness is defined by 


v (Oi - Ci ) 2 


*“d n-p- 1^' cr 2 

r 1= 1 l 


( 1 ) 


where O, is the relative flux measured at the wavelength A,-, C , 
is the relative flux computed from the model at the same wave¬ 
length Ai, crj is the variance of the data point i that is defined as 
the (S/N) 1 ratio at the same wavelength An is the number of 
data points, and p is the number of free parameters of the model 
so that n — p — 1 is the number of degrees of freedom (d.o.f). 

For the spectra of HD 10700 and HD 122563, we did not 
have the cr, at the wavelength A,. Although it is known that the 
S/N depends on the number of recorded photons I (~ V7), we 
adopted the mean constant S/N ratio for these stars because we 
only fit the wings of each line up to 0.9 where the S/N ratio 
is not changed dramatically. For Procyon, for example, the rel¬ 
ative flux at a S/N of 0.85 is 10% lower than that with a S/N 
of 0.995. Apart from the inner wings, which have less strongly 
blended lines because they are saturated by the hydrogen line, 
equal weights were assigned for all spectral windows of a line 
profile. 

There were two free parameters, otcma and the effective tem¬ 
perature, while the other parameters (logg, [Fe/H], /') were fixed. 
Assuming the Poissonian error distribution, which approaches 
the normal distribution at high S/N, the minimization of x 2 pro¬ 
vides the maximum-likelihood estimate of the parameters, in our 
case for qtcma anc * Teg. 

We also compared a grid of Balmer-line profiles for a given 
model grid with a 25 K interval in T e g- and a 0.1 interval in o'cma- 
By varying the free parameters of the model, we computed the 
value of x 2 at each step. First we fixed one of the two free pa¬ 
rameters, T e ff, and analyzed the function of x 2 , which smoothly 
changes with varying o'cma and follows an asymmetric parabolic 
form with one single minimum. The best otcma and 7/ f | were 
chosen from the requirement of the lowest x 2 value. The vari¬ 
ation of temperature within the error bars allows us to estimate 
the upper and lower values for our o'cma- An example of the fit¬ 
ting for H ff profile in the metal poor star HD 103095 is shown in 

FigOU] 


6.1. Errors 

Estimated errors are presented in Table [2] for the program 
stars. The strength of the hydrogen line depends only weakly 
on the surface gravity or the metallicity. We considered how 
the typical errors (0.1 dex) for either of them may affect the 
result. Microturbulence does not affect the formation of the 
Balmer-line wings (IFuhrmann et al.|[l993l) . therefore we did not 
consider the influence of micro turbulence. Only the tempera¬ 
ture can strongly affect the wings of hydrogen lines. With in¬ 
creasing temperature the wings of hydrogen lines tend to be 
stronger. The effective temperature is the main source of er¬ 
rors in our results. The influence of other sources of error, such 



Fig. 10. Example of fitting the H ff wings in the metal-poor star 
HD 103095. Bold dots correspond to the observed data. The theoreti¬ 
cal H„ profile was calculated with (dashed curve) and without (solid 
curve) blending lines. Only bold dots were used for determining the^ 2 
statistic. 


as Stark-broad ening, self-broa d ening , and H e-br oadening, were 
consid ered in iFuhrmann et al.1 (1 19931 |T994) and Barkle m et al.1 
(2002j). The latter showed that the influence of Stark- and self¬ 
broadening on the behavior of Balmer lines is insignificant for 
solar metallicity stars, but tends to be stronger for metal-poor at¬ 
mospheres. Moreover, the determination of the continuum place¬ 
ment is one of the most significant uncertainties in this work. 
For the sample spectra, the error in the continuum placement 
of the high-quality spectra is approximately 0.5% of the contin¬ 
uum flux, depending on the S/N ratio. As a result of the inter¬ 
dependency of errors, it is difficult to precisely account for the 
combined effect, therefore we divided all errors into two groups: 
those propagated from the stellar parameters (7^ c n, log g, and 
[Fe/H]), and those caused by observational aspects (continuum 
placement and fitting error). 


7. Results 

In Fig. QT]we plot the profiles of the Balmer lines for three of 
our program stars (HD 122563, HD 103095, and the Sun). The 
profiles are only slightly different under LTE or NLTE assump¬ 
tions for dwarfs, but exhibit significant variations with a cma on 
their wings. The changes with cccma from 0.5 to 2.5 are simi¬ 
lar (for H„) and even larger (for H /; and H y ) to that of T e ff of a 
few hundreds of Kelvin. For the metal-poor giant HD 122563, 
the NLTE effect plays a more important role than o'cma for H<> 
profile, but is less sensitive than o'cma parameter. Figure QT| also 
indicates that if the q:cma parameter in the model atmosphere is 
not correct, the T e g may be biased by up to a few hundreds of 
Kelvin by fitting Hg or H r line profiles. Therefore, caution must 
be exercised when choosing ctcma in ID model atmospheres. 


7. 1. The Sun 


The Sun is the best-known star with high-quality observations 
and accurate physical parameters. Following iFuhrmann et al.l 
(tl99l . we attempted to calibrate qtcma for the so l ar atm o- 
sphere. Our work differs from that of IFuhrmann et al. ( 19931) in 
that we have used the self-broadening theory of Barklem et al.l 
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Table 2. Estimated errors in CMA factor for H„ profiles. 


Error 

Sun 

Procyon 

HD 10700 

HD 103095 

HD 39587 

HD 6582 

HD 217014 

HD 22049 

HD 122563 

HD 45282 

HD 140283 

ATeff 

±0.2 

±0.4 

±0.3 

±0.3 

±0.4 

±0.3 

±0.3 

±0.3 

±0.1 

±0.3 

±0.2 

A log g = 0.1 dex 

- 

- 

- 

±0.1 

- 

±0.1 

- 

- 

±0.1 

±0.1 

±0.1 

A[Fe/H] = 0.1 dex 

- 

±0.2 

±0.1 

±0.1 

- 

±0.1 

±0.1 

±0.1 

±0.1 

- 

- 

Fitting procedure 

±0.2 

±0.2 

±0.1 

±0.2 

±0.1 

±0.5 

±0.1 

±0.1 

±0.2 

±0.2 

±0.2 

Continuum 












(0.5%) 

±0.3 

±0.5 

±0.2 

±0.3 

±0.1 

±0.3 

±0.1 

±0.1 

±0.2 

±0.2 

±1.0 

Total 

±0.4 

±0.7 

±0.4 

±0.5 

±0.4 

±0.7 

±0.4 

±0.4 

±0.4 

±0.4 

±1.1 


HD122563 



AA (A) 


HD103095 sun 



AA (A) AA (A) 


Fig. 11. Comparison of the line profiles of H a , H^, and H y (from bottom to top) for HD 122563, HD 103095, and the Sun. The x-axis is the 
wavelength shift relative to the center of the Balmer lines. The colors are coded with different « C ma (from 0.5 to 1.5), while the solid lines and 
dashed lines represent profiles using LTE and NLTE line formations, respectively. The dotted lines are profiles with LTE models of a cma = 0.5 
but T cff + 100 K. For clarity, offsets of 0.2 and 0.4 are added to the Hg and H y , respectively. 


(12000}). For investigations of the Sun we employed the solar at- 
las ( iKurucz et al.lfl984t) . We adopted an effective temperature of 
5777 K with an error bar of +40 K. By fitting theoretical H„ pro¬ 
files to reduced observational profiles, the two parameters T e g 
and (+cma were varied within error bars and reasonable physical 
limits, respectively. The results are shown in Fig. [12] where the 
X 2 minimum indicates the best T e g and qtcma- 

Figure IT3l compares the left wing of H„ line profile for the 
Sun with the theoretical spectrum for ocma = 2.0 (solid line). 


The combination of T e ff=5770+20 K and i+cma = 2.0 ± 0.2 is 

preferred based on the degree of agreement between the obser¬ 
vational and theoretical H„ profiles. Nevertheless, it is important 
to note that the value of x~ is still high, which suggests that our 
agreement between the observational and the theoretical H„ pro¬ 
files is not perfect. This means that the quality of the observed 
solar flux spectra is higher than the quality of theoretical surveys. 

From comparison with SED, a low value for ctcma < 10 
is excluded. The lack of fitting in the region below 5000 A does 
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not allow precisely determining the convective efficiency, but the 
measurements are consistent with an ctcma of 2.0. 


7.2. Procyon 


Procyon (HD 61421, HR 2943) is one of the brightest stars in the 
night sky (V = 0.34). It is a binary system consisting of an F5IV- 
V star (Procyon A) with a mass of 1.497 ± 0.037 M e and a faint 
white dwarf (ProcyonB) with a mass of 0.602_± 0.015 M 0 , or¬ 
biting each other with a period of 40.8 years (iGirard et al.i2000i) . 
Spectroscopic analyses gave rather discrepant effective tempera¬ 
tures, from 6470 K (by fitting the Baimer lines; Zha o & Gehrenl 
120001) to 6850 K (by measuring the equivalent widths of a 


large set of iron lines; iHeiter & Luckl 20031). The prox imity 


of Procyon (r = 284.56 ± 1.26 mas; Ivan Leeuwe 
the Earth allows us to directly measure its 


; proxtm 
d l2007i) 


to 


ter (e.g., Hanburv 

Brown et all 1974; Mozurkewich et alJ 1991; 

Aufdenberg et al.1 

2005; Chiavassa et al. 2012). Kervella et al. 


agrees well with othe r works (e.g.. lAllende Prieto et alll2002t 
Fuhrman n et aljfl997l> . We used a high-quality spectrum of H a 
for Procyon and applied the same fitting method as described in 
the previous sections. 


The results are displayed in Fig. [T2l where the y 1 minimum 
indicates the best T e g and ctcma- For this star we find 77 ir = 
6524 + 50 K and otcma = 1.9 + 0.7. The error of otcma is not 
symmetric in this case, but we adopted the largest. The left wing 
of the H,, line profile for Procyon is shown in Fig. [13] where two 
theoretical spectra with ctcma = 1.5 (solid line) and otcma = 2.3 
(dashed line) are plotted as well. 

The value preferred from the SED fitting for »cma ~ 1.4, 
nevertheless, there are large uncertainties with respect to opac¬ 
ity data. As a result of the uncertain UV flux comparison, the 
values given by SED fitting and that by Balmer profiles do not 
contradict each other. 


7.3. Program stars 

The final results of our fitting are given in Table Q] including 
the names of the program stars, their parameters, the best o-cma 
along with the limits, and the^ 2 values that indicate the best fit in 
each case. We also present the best-fit T e g- within the error bars 
that was obtained during the x 2 minimization procedure. Tem¬ 
peratures with asterisks are the best-fit temperatures. The results 
from Hg and H y for some stars are less reliable. Even though 
good-quality observations of Hg and H y are available for all re¬ 
maining stars, these lines were not adopted because choosing a 
reliable mask was problematic. 

The observed spectra and the best-fitting line profiles from 
this method are plotted in Fig. [14] When applying the above- 
mentioned method, we have to keep in mind that weak blend¬ 
ing lines may, nevertheless, be presented in a chosen windows. 
Therefore, it would be better to inspect the best profiles visually 
to give preference to the case where all observed points are be¬ 
low the theoretical line profile, but the x 2 ~ method minimizes 
the deviations without considering whether the observed points 
are above or below the theoretical profile. In this context, for 
the x 2 method, »cma will be slightly underestimated, while the 
best-fit T e ff will be overestimated. 

The derived T eS for stars HD 10700, HD 39587, HD 217014, 
HD 22049, HD 122563, HD 45282, and HD 140283 agree well 
with those derived from astrometric and IRFM methods. How¬ 
ever, for HD 103095 we were unable to find an appropriate 77 ir 




Fig. 12. Survey of the reduced x 2 obtained from profile fitting for 
the observed solar flux spectra (top panel) and spectra of Procyon (bot¬ 
tom panel) with different Ucma and TefS- The preferred parameter area 
is located where x~ is minimal. All y 2 are not normalized and indicate 
a good fit in each case. 


in the temperature range T e g = 4820 ± 100 K. Thus, for this star 
we adopted a higher spectroscopic temperature of 77 ir = 5110 
K. 

The individual results for our stellar sample are reported be¬ 
low. 

HD6582 (//Cas, 30Cas, HR321) is a close binary sys¬ 
tem that consists of a slightly metal-deficient subdwarf (com¬ 
ponent A, V = 5.17) and a component B fainter by 5.3 mag¬ 
nitudes (V = 10.5), with a separation of only 0".40~1".40 
dMason et al.ll200l]) . The spectroscopic parameters adopted here 
for HD 6582 A (T eS = 5387 K, logg = 4.4 5. fFe/H] = -0.83 , 
and £ = 0.89 kms ~*) were taken from Fuhrmann (2004). 
Bov aiian et alJd2012l) measured its limb-darkened angular diam¬ 
eter Old = 0.972 ± 0.009 mas with CHARA, corresponding to 
the stellar radius of 0.790 ± 0.009 R e , combining with the Hip- 
parcos parallax (n = 132.38 ± 0.82 mas: Ivan LeeuwenH2007i) . 
This is agrees well with the spectroscopic radius ( R = 0.8 R e ), 
but the interferometric 77 n = 5264 K is ~120 K lower than that 
of lFuhrmann i (2004). By comparing the H, ( profile, we found the 
best o-cma = 2.3 for spectroscopic parameters, while acMA =1.7 

for interferometric parameters. Despite this large discrepancy, 
the fitting method with the interferometric model (b) led to lower 
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Fig. 13. Left wings of LL line profiles for the Sun with o-cma = 2.0 (top 
panel) and Procyon (bottom panel), where two theoretical spectra with 
Ocma = 1.5 (solid line) and ocma = 2.3 (dash line) are shown. 


X 2 and supports a higher 7^ c fi = 5370 K. This result supports 
the relative high 7Tir of 5339 K derived by the photom etric 
b — y = 0.437 and ci = 0.213 (Hauck &Mermillio(li 19981) with 
the color-7 e ff relation of lAlonso et alJ (119961) . 

HD 10700 (rCet, 52Cet, HR 509) is a nearby and inac¬ 
tive G8 V star, with a distance of only 3.650 ± 0.002pc {n = 
273.96 ± 0.17 mas: Ivan Leeuwetil 12007 ). It is one of the most 
frequently observed targets in several precise radial velocity 
surveys for extrasolar planets, and iTuomi et al.l (1201 3b reported 
five planets with minimum masses of up to 6.6 M® orbiting it. 
Spec troscopic analyse s indicate that its 7 W ranges from 5283 
K (1 Valenti & Fischeil 120051) to 5420 K dTakeda et alJ 120051) . 
and [Fe/H] range fr om -0.43 dTakeda et al. 2005) to -0.58 
dLuck & Heitejl2005l) . Interferometr ic measurements with VLTI 
gave Old = 2.078 ± 0.031 mas dDi Folco et alJ 120041) . and 
7 e ff was perfectly consistent with the spectroscopic parameters 
(7 eff = 5377 K, logg = 4.53, [Fe/Hl = -0.49, and 7 = 0.8 
kins -1 ; iMashonkina et al.l 1201 ll) . With a least x 2 fitting, we 
found the best a cma = 1-8, and the best T e s = 5330 K agrees 
with the T R ff — 5346 K inferred from the b - y and c \ color 
dHauck & Mermilliodlll 998t lAlonso et al.l[l996h . 

HD 22 049 (eEri, 18Eri, HR 1084) is a nearby (n = 310.94 ± 
0.16 mas; Ivan Leeuwenl 2007), c hromospherical ly active K2 V 
star orbited by a dusty ring ( Greaves et al.ll 1998b . Atmospheric 
parameter determinations based on high-resolution spectra gave 
very consistent results, w ith a T e ff ranging f rom 5054 K 
(lFuhrmannl2004 )_to 5200 K dLuck & Heite 1^20051), log g r angin g 


from 4.40 dMishenina et al.ll2013l) to 4.7 2 (flhkeda et al.ll2005b . 

and [Fe/Hl ranging from -0.18 dGhezzi et al.l 120 1 oh to +0.06 


(Take da et alJl2005h . We adopted the parameters of[Fuhrmann 
(2004!) (T„ ff = 5054 K, log g = 4.5 4, [Fe/H] = -0.09, and 7 = 
0.9 kms' 1 ). iDi Folco et al.l ( 2004 ) suggested an interferomet¬ 
ric angular diameter of Old = 2.148 ± 0.029 mas and an ef¬ 
fective temperature of T e ff = 5122 K. H,, profile fitting with 
spectroscopic and interferometric methods result in quite con¬ 
sistent o-cma = 2.2. The best-fit T e ff = 5045 K also agrees 
with the photometric T P ff = 5074 K from the b - y and c\ 
dHauck & Mermilliodfilim 

HD 39587 (x 1 Ori, 540ri, HR 2047) is a variable star of RS 
CVn type, orbiting by a faint (Va = 4.4, Vb = 8.8) companion 


with a mass of 0.15 + 0.01 M® dlrwin et al.lll992l) . lBovaiian et al.l 
d2012l) measured its angular diameter Old = 1-051 + 0.009 mas 
with a precision of 0.9%. This yield a radius of 0.9 79 + 0.009 R n 
with a Hipparcos parallax o f 115.43 ± 0.27 mas dvan Leeuwenl 
2007). Boyaiian et al. (2012) also gave its effective temperature 
as 7 e ff = 5961 K, which generally a grees with the spectroscopic 
temperature reported bv iFuhrmannl d2004l) . We found the best-fit 
<+cma =1.5 for spectroscopic parameters and a cma = 1.9 for 
interferometric parameters. 

HD 122563 (HR 5270) is the most metal-poor halo star 
([Fe/H]~ -2.6) that can be seen by the naked eye (V = 6.2) 
in the night sky. It is the only giant i n our program stars and 
therefore provides unique constraints. iMashonkina et al.l (1201 ll) 
presented a detailed analysis of the NLTE formation of iron lines 
for this star, and we used the same parameters (7 e ff = 4600 K, 
logg = 1.60, [Fe/H] = -2.60, and £ = 1.95 kms -1 ). Recent inter¬ 
ferometric measurements with CHARA and PTI dCreevev et al.l 
12012|) gave its angular diameter as Old = 0.948 ± 0.012 mas, and 
Left = 4598 ± 41 K, which is consistent with its spectroscopic 
temperature. We found the best-fit <+cma = 1.0 with H„ profile, 
while the fitting of H/- ; and H y gave satisfactory agreement with 
<+cma = 0.5. The + cm a is l° wer than that of the dwarfs, implying 
that the stellar convective energy transport become less efficient 
after stars evolved off the main sequence. 

H D 217014 (51 Peg, H R 8729) is a nearby (n = 64.07 ± 0.38 
mass; Ivan Leeuwenl 120071) G5 V star with t he first extrasolar 
plane t ever discovered orbiting a Sun-like host (Mavo r"& Ouelozl 
119951) . Spectroscopic anal ysis gave its effective t emperatures 
as ranging from 571 0 K dMaldonado et al.l 1201 2l) to 5832 K 
dRamtrez et al.ll2009h . The narrow- band photometric color in- 
dices b — y — 0.416 and c\ = 0.371 dHauck & Mermilliod ,1998) 
indicate a rather low T e g - 5661 K ( Alonso et al.lll996 ). The 
parameters adopted here (T . ff = 5793 K, log g = 4.33, [Fe/H] = 
+0.20, and £ = 0.95 km s~ 1 . IFuhrmannl 19980 agree for the effec¬ 
tive temperature with the interferometric results^of T e ff = 5804 
K and Old = 0.748 ± 0.027 mas dBaines et al.ll2008l) . We found 
a best-fit of ocma = 2.0 with both parameters giving consistent 
results. 

HD45282 is a GO metal-poor star ([Fe/H] = -1.5) at low 
Galactic latitude (b = -3.9°). Previou s spe ctroscopic analyses 
gave its T e s as ranging from 5150 K (iFulbrightil2000l) to 5344 K 
dGratton et al.ll2000h . Consid ering its high reddeni ng with E(B - 
V ) =; 0.82 at thi s direction ( Schlegel et al. 1998), we used the 
IRFM T e ff from ICasaerande et al.l d2010h . which is thought to 
be less affected by interstellar medium. The Balmer profile x 2 
fitting gives i+cma =0.1 for H,, and ctcma = 0.5 for H /; with 
the best-fit T e ff = 5310 K. The error of 0.05 on E(B - V) causes 
an uncertainty of ~ 300 K on T e ff, therefore the photometric 
temperature is not reliable. 

For the Sun we used the widely adopted parameters of T e ff = 
5777 K, logg = 4.44, [Fe/H] = 0.0, and /■' = 0.90 kms- 1 and 
for Vega T eff = 9550 K, logg ^ 3.95 , [Fe/H] = -0.50, and = 
2.00 kins -1 dCastelli & Kuruczll 19941) . The solar o-cma from the 
H,, Balmer line is 2.0, while a determination based on the solar 
flux suggests a value of 1.5-1.75. Although this is consistent 
within the error bars, it has to be noted that, for both the Sun 
and Procyon, SED fitting leads to a notably lower value of <+cma 
than Balmer-line fitting. 


7.4. o-cma 

Figure[15] shows our final results for the efficiency parameter 
<+cma as determined from Balmer-line fitting. Except for evolved 
stars, acMA is consistent with a value of about 2.0 for all objects 
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Fig. 15. Determined Ocma f° r the program stars. An error of q-cma of ~ Fig. 16. Difference between our determined ucma and the one predicted 
±0.4... 0.7 maybe assumed for all stars except HD140283 (T e ff = 5780 bv lMaeic et al.l J201 4| > . 

K, logg = 3.70 where the error reaches » ±1.1. 


within the given insecurities. There is no obvious correlation 
with metallicity and a weak trend of increasing a for the coolest 
main-sequence objects in the investigated range of F- and G-type 
stars. This confirms the behavior and numerical values predicted 
bv lMagic et alJ (1201 4l) for this part of the HRD. 

For the evolved stars of our sample otcma is of about 1.0. 
Although iMagic et alJ ( 2014 ) predicted the same trend of lower 
convective efficiency for objects with lower log#, the magni¬ 
tude of the change is much greater in our sample than that sug¬ 
gested by the 3D simulation of lMagic et al. ( 2014 ). Our sample 
mainly consists of main-sequence stars and only contains three 
evolved objects with well-determined parameters. Nevertheless, 
the stars clear I v_s how a behavior different from that suggested by 
IMagic et~ahl (120141) . The three objects span a wide range of tem¬ 
peratures, therefore the evolutionary state seems to be the main 
discriminator that separates them from the rest of the sample. 


8. Discussion and implication 


Figure[l6] shows the difference between the suggested ctcma and 
that from IMagic et al.l ( 2014 ). who adopted a calibration using 
the entropy jump (Fig. 3 in their work). As described above, the 
differences are small along the main sequence and increase dra¬ 
matically later in the evolutionary sequence, which suggests that 
there is either a stronger variation of »cma as predicted by the 
STAGGER grid, or that our modeling of the Balmer lines is in¬ 
correct. The latter is excluded because we generally determined 
a temperature that agreed well with astrometric and IRFM data. 

Based on the similar method of FAG93, our analysis indi¬ 
cates that if the theory of Ba rklem et alJ ( 2000) for Bal mer-lin e 
broadening is adopted instead of that of lVidal et al.l (11 97(1 11973l), 
while the old ODF version is replaced by the MAFAGS-OS ver¬ 
sion, it is very unlikely to derive a low value of »cma = 0.82. 
For F- and G-type main-sequence stars, qtcma ~ 2.0 is derived, 
and the value decreases toward 1.0 for more evolved stars. 


A variation of the convective efficiency naturally influences 
the temperatures determined by the Balmer-line method. It fur¬ 
thermore changes stellar SEDs and in turn filter fluxes and color 
determinations. All stellar lines that have contribution functions 


down to the convection scone (very strong lines) are subject to 
potential changes as well. 
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Fig. 14. Best fits of the lines in our stellar sample. Solid lines are the theoretical spectra corresponding to the best fit and dotted lines represent the 
observed spectra that were used for determining the^ 2 statistic. The horizontal dotted lines show the limit of reasonableness of the approximation 
and allows seeing residuals. 
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